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r T , 1 Abstract 



Effects of non-resonant photoproductions arising from two different ttNN couplings are investi- 
gated in the — ► ttttN reaction. We find that the pseudoscalar (PS) ttNN coupling is generally 
\ preferable to the pseudovector (PV) irNN coupling and particularly the total cross sections are 

o ' 

successfully described by the model with the PS ttNN coupling. In order to see the difference be- 
■ tween the two couplings, we also show the results of invariant mass spectra and helicity-dependent 

O 

cross sections in various isospin channels calculated with the PS and PV couplings. 
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I. INTRODUCTION 



Many experiments [l-9]of two-pion photoproductions on the nucleon have been performed 
in the past over the second nucleon-resonance energy where the iV*(1520) resonance plays 
an important role. Two charged pion photoproduction, i.e., yp — > n + 7i^p, is well understood 
by theoretical models. However, none of the theoretical models have succeeded to explain 
the data of all isospin channels simultaneously. This means that there are some unknown 
production mechanisms which are not taken into account yet. 

Recently, the total cross sections and invariant mass distributions of two-pion photo- 
productions followed by neutral pion emission, i.e., yp — > 7r + 7r°n[6], yp — > 7r°7r°p[8] and 
7n — > 7r~7r°p[9] reactions, have been measured at the Mainz accelerator facility MAMI using 
the detector system with improved resolution at the photon energy up to around 0.8 GeV. 
The aim of these experiments was to obtain the information on the structure of the nucleon 
resonance and explore the reaction mechanisms. The •yp — > n + 7i°n and yn — > 7i^n p re- 
actions have attracted special attention, since the detailed study of them could provide a 
new aspect on the reaction mechanisms which are related to p meson production. There 
are notable differences between the ir + ir° (it~it ) and 7r°7r° photoproductions. The p meson 
production as an intermediate process is allowed to the 7r + 7r°( tv~ti ) photoproduction but 
is forbidden to the 7r°7r° photoproduction due to isospin conservation. On the other hand, 
the isospin I = mv system such as er meson may contribute to only the 7r°7r° photopro- 
duction. Furthermore, the strength of the A Kroll-Rudermann process in the 7r + 7r°(7r~7r°) 
photoproduction is weak compared with that in the 7r + 7r~ photoproduction where its pro- 
cess dominates and is suppressed in the 7r°7r° photoproduction. Based on these characteristic 
features and the comparison between the measured n + n and 7r°7r° invariant mass distribu- 
tions, W.Langgartner et al. [6]concluded that the p decay of the iV*(1520) resonance was 
directly observed in the yp — > 7r + 7r°n reaction. 

Theoretically, several works [10-14] have been already done to explain the total cross 
sections of two-pion photoproductions in various isospin channels. The total cross sections 
of the 7p — > -K + n~p reaction have been well reproduced by several theoretical models[10-14] 
where the A Kroll-Rudermann, A-pion-pole and iV*(1520) resonant terms were found to be 
dominant. These models, however, could not predict the total cross sections of the above 
reactions accompanied by neutral pions consistently. The models of Tejedor-Oset[10] and 
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Murphy-Laget[ll] underestimated the cross sections of yp — > n + n°n seriously, although the 
former model could reproduce the cross sections of yp — > n -n°p. 

The model of Ochi-Hirata-Takaki[12, 13] could reproduce the yp — > 7r + 7r°n and yn — > 
7r _ 7r°p cross sections as well as yp — > it + it~p cross sections, although it underestimated 
the 7p — > 7i°n°p cross sections. Their calculations indicated that the p meson productions, 
i.e., the iV*(1520) — > pN process and p Kroll-Ruderman process, play an essential role 
in the •yp — > n + n°n and yn — > 7r~7r°p reactions where the mass of p meson produced in 
the intermediate state is always smaller than the on-shell value at the relevant energies. 
In their model, the p meson is treated in a dynamical way where the finite-range pirn 
form factor is assumed. In order to obtain large cross sections of the yp — > 7i + 7i°n and 
7n — > 7T~7T°p reactions, a rather soft pnir form factor, which makes the contribution of the p 
Kroll-Ruderman term large, was needed. Because of the soft pirn form factor, however, the 
dynamical model for the p meson overestimates the nn p-wave (isospin I — 1) phase shifts at 
low energies. They speculated that the large p Kroll-Ruderman term in their model might 
simulate a background process in the isospin 1 = 1 channel effectively rather than the p 
meson production. The presence of such a background process is inferred from the fact that 
in the isovector 7nr spectral function derived from the 7iir — > NN helicity amplitudes[15], 
there is a strong enhancement near the nn threshold as well as the resonant structure by the 
p meson. This bump at low energies is actually due to the non-resonant process described 
by the partial amplitude of the nucleon Born term projected to the I = J = 1 irir(NN) 
channel. 

In this paper, motivated by the above speculation, we will investigate the effect of the 
non-resonant reaction mechanisms, particularly, the background terms arising from the 7rNN 
coupling that were so far considered to have only small contributions to the cross sections 
and were always taken to be the pseudovector (PV) coupling. In the studies of single 
pion photoproductions, it has been shown that the PV 7riViV coupling is preferred at low 
energies but the pseudoscalar (PS) reNN coupling is needed to get a better description at 
higher energies[16] and furthermore the two irNN couplings lead to rather different cross 
sections for the neutral pion production [17]. These facts imply that the PS coupling becomes 
important at larger off-shell nucleon momenta and thus will have a significant influence on 
the two-pion photoproductions accompanied by the neutral pion such as the •yp — > 7T + n°n 
and 7p — * 7r°7r°p reactions. In our calculations with two different couplings, we will show 
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that the non-resonant photoproduction by the PS coupling significantly contributes to the 
total cross sections of the two-pion photoproductions involving the neutral pion compared 
with that by the PV coupling and consequently plays a similar role with the strong p Kroll- 
Ruderman term introduced in the model by Ochi et al. [12, 13]. From the comparison between 
the full calculations with resonant processes and the experimental data, furthermore, we will 
demonstrate that the PS coupling is more favored than the PV coupling for the two-pion 
photoproductions at the relevant energies. 

Recently, Nacher et al.[14] have improved the model of Tejedor-Oset[10] by including the 
A(1700) production and the p meson effect arising from the iV*(1520) production. In their 
calculations with the PV coupling, they found that the p meson effect largely increased the 
7P — > ii + 7i n total cross sections compared with the previous model and put the calculations 
close to the data, although there still remained some disagreement with the data around 
the peak for both the jp — > 7i + n n and ^p — > n 7i p cross sections. We will discuss why 
their model can largely improve the calculations of the •yp — > 7i + 7i n cross sections without 
introducing the strong p Kr oil- Ruder man term used in the model by Ochi et al.[12, 13]. 

In section II, we will discuss the background processes and ttNN coupling. In section 
III, we will review how to treat the resonance processes. In section IV, we will show our full 
calculations of the total cross sections, invariant mass spectra and helicity-dependent cross 
sections and discuss the difference between the PS and PV couplings from the comparison 
with the data. In section V, we will give our concluding remarks. 



II. NON-RESONANT PROCESSES AND vriViV COUPLING 



In this section we discuss the non-resonant processes arising from the ttNN coupling and 
vector mesons. There are two types of nNN couplings: the pseudoscalar (PS) coupling and 
the pseudovector (PV) coupling. The Lagrangian is 



£*NN = WkNN i>l b T lf)4>, (1) 

for the PS coupling, and 

C% N = ^^7 5 r W^, (2) 
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for the PV coupling, where ip and <f> are the nucleon and pion fields, respectively, gl NN /47r = 
14.4 and fnNN/mn = —QttNn/^Mn and m n and Mn are the pion and nucleon masses, 
respectively. These couplings are equivalent when both the nucleons are on-shell. If a 
photon line is attached to the 7iNN system so as to become gauge-invariant, one obtains the 
Feynman diagrams of the Born terms for the one-pion photoproduction whose expressions 
depend on the ttNN coupling. The 7r Kroll-Ruderman term is included in the PV Born 
terms. The Born terms calculated with the two couplings are rather different each other 
because of their different off-shellness. The PV coupling is preferable to the PS coupling 
at low energies because the PV Born terms are consistent with low energy theorems and 
current algebra predictions. In fact, the model with the PV coupling is able to reproduce 
the E 0+ multipole up to the A(1232) resonance energy region but the model with the PS 
coupling fails. As the incident photon energy increases over 500 MeV, however, the pure 
PV coupling cannot explain the E 0+ and M x ~ multipoles and the PS coupling is needed to 
describe them[16]. This suggests that the irNN vertex for the far off-shell nucleon is largely 
pseudoscalar in nature. We note that only the multipoles £"0+ and Mi_ are affected by 
changing the coupling scheme and the two ttNN couplings give rise to significantly different 
cross sections for the neutral pion photoproduction but are almost indistinguishable in the 
charged pion photoproduction. Therefore, it is interesting to see the difference between 
the PV and PS couplings in the two-pion photoproduction around the iV*(i520) resonance 
energy region where the far off-shell nucleons are involved in the intermediate state. 

Before going to the — > n-nN reaction, we discuss briefly the strong interaction part 
relevant to the two-pion production. In the nN — > nN scattering at low energies, the Born 
terms constructed with the PV coupling are also preferable to those with the PS coupling 
like the 7V — * nN reaction. In this case, for instance, the s-wave isoscalar scattering length 
calculated from the PS Born terms is very large and in disagreement with the data. In 
dispersion relation theory, the PS Born terms correspond to nucleon-pole terms and the PV 
Born terms are understood to include strong corrections coming from the dispersive integral 
in addition to the pole terms. These corrections are important to well describe the low 
energy data of the nN scattering and may be partially related to the a meson exchange 
contribution as inferred from the linear a model. When constructing the model for the 
•yN — > ixixN reaction with the PS nNN coupling, these corrections must be taken into 
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FIG. 1: Diagrams of Born terms arising from itNN coupling ((a)-(g)) and the effective Lagrangian 
ofEq.(3) ((h)-(k)). 

account. To do so, we introduce the following effective Lagrangian 

C^ NN = ff^</> ■ <f> + ^§^^ T x d„d>). (3) 

The sum of the PS Born terms and their correction terms calculated from the above effective 
Lagrangian are equivalent to the PV Born terms for the strong interaction processes such 
as 7iN — > 7iN and N — > ttttN. We note that the 7rN isoscalar scattering length calculated 
in these models are consistent with the value obtained by taking m a — > oo in the linear a 
model, where m a is the a meson mass. 

The diagrams shown in Fig.l for the 7iV — > 7r7riV reaction are obtained by attaching an 
external photon line to the diagrams of the N — > tiitN processes based on the requirement 
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of the gauge invariance. The diagrams calculated with the PV coupling correspond to 
Fig.l(a)-(g) and especially the diagrams (f) and (g) include the ir Kroll-Rutherman term 
arising from the derivative nNN coupling. On the other hand, the diagrams calculated 
with the PS coupling and the above effective Lagrangian correspond to Fig.l(a)-(e) and 
(h)-(k), respectively. They will be referred to as PV model and PS model, respectively. For 
the second term (isovector term) of Eq.(3), only the diagram (k) of Fig.l is computed in 
actual calculations because the contributions of the other diagrams are negligible. Here the 
Lagrangians for the 7iVW and 77r7r couplings are given by 

C lNN = -e^FrfAp - F 2 a^(d„A^, (4) 
£ 7 ™ = e(d M </> x (f)) 3 A^, (5) 

where A p is the photon field, and F x and F 2 are the electromagnetic form factors which are 
taken to be F x = 1 and F 2 = 1.79/(2M N ) for the proton and F x = and F 2 = -1.91/ (2M N ) 
for the neutron, respectively. The other Lagrangians for coupling of a photon and hadrons 
are obtained from the minimal substitution (<9 M — ■> + ieA^) in the Lagrangians of 7iNN 
and nnNN derivative vertices. The resultant amplitudes for the — > irirN reaction will 
give rather different cross sections depending on whether the PS model or PV model is used 
because of their different off-shell behavior, as expected from the above discussion about the 
7iV — > irN reaction. We will show the difference from the numerical results later. 

As other non-resonant processes, we include the contributions of vector mesons such as 
p meson and uj meson. In the analysis of the 7V — > nN reaction, these contributions 
are known to be important although their magnitude is small. The hadronic Lagrangians 



involving these vector mesons are 

W = - i= (g v pNN r - f|^^ p» ■ r ^ (6) 

^ = -/^-(0x^), (7) 

£u,nn = - ? (gl NN Y - g ikr^ v d^ ^ (8) 



where and uo p are the p and u meson fields, respectively and the coupling constants 
are taken to be g^ NN = 2.9, gJ NN = 18.15, f p = 6.0, g% NN = 7.98 and g^ NN = 0. The 
electromagnetic Lagrangian for the 77r u; coupling is 
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FIG. 2: Diagrams of non-resonant processes arising from p and 10 mesons. 



= - ^ w (^)7r d P K), (9) 
where g 17TW = 0.374 e and = 783MeV. The other electromagnetic Lagrangians are 
derived from the minimal substitution as done for the PS and PV models before. The 
diagrams for the — > txtiN reaction involving the p and u mesons are shown in Fig. 2 and 
are calculated by using the above Lagrangians. For the u meson contribution, diagrams (g) 
and (h) in Fig. 2 are taken into account. The diagram (c) (called p Kroll-Ruderman term) 
comes from the derivative pNN tensor coupling. In the diagrams of (a)-(d) in Fig.2, the 
p meson decays into two real pions directly and so the p propagator D p must include the 
decay effect, whose form is assumed to be 

^p(v^) = : -5— L-F7 7^v ( 1Q ) 



m 



with 

2/, 2 



UVs) = 5^qL (ID 



Here the p meson mass m p = 775 MeV, its width Y p (m p ) = 150 MeV, and q cm stands 
for the pion momentum in the rnr center of mass system. In our present model, the p 
meson displayed in Fig. 2 is treated in the same way used in Ref.[10, 14] where the pixix 
vertex function is simply given by the Lagrangian of Eq.(7). Even though one uses the pmr 
vertex function with the finite-range form factor employed in Ref.[12, 13], there is no drastic 
numerical change from the above way in the magnitude of the p Kroll-Ruderman term as 
far as one uses the range parameter determined so as to reproduce the nn p-wave (isospin 
1 = 1) phase shifts. 

For the diagrams involving the off-shell meson coupled to nucleon in Fig.l and Fig. 2, we 

take into account a form factor of the following form: 

. 9n A 2 — m 2 . . 

F{q) = A^-T' (12) 

where m and q 2 are the meson mass and the square of four-momentum, respectively and 

the range parameter A is taken to be A = 1.25 GeV for ir, A = 1.4 GeV for both p and uj, 

respectively [14]. The pion form factor is used for the diagrams (d) and (e) in Fig.l and the 

p meson and u meson form factors are used for (a)-(f) and (g)-(h) in Fig. 2, respectively. For 

the diagrams (f) and (g) in Fig.l (rr Kroll-Ruderman term), the same pion form factor is 

used and evaluated at the momentum transfer between the incident photon and the outgoing 

pion. The diagram (c) in Fig. 2 (p Kroll-Ruderman term) is treated in the same way as the 

7i Kroll-Ruderman term. We note that the gauge invariance of the transition amplitudes 

is destroyed by the inclusion of the decay width in the propagator and the form factor for 

the hadronic vertex. However, we consider such strong interaction corrections are more 

important than the gauge invariance from a phenomenological point of view. In fact, the 

form factor used in the A Kroll-Ruderman term (see Fig.5(a)) influences significantly on 

the magnitude of the 7p — > 7r + 7r~p cross section and provides a good agreement with the 

data[10]. 

Now we discuss the differences between the PS and PV models through numerical cal- 
culations of the cross sections. The cross section for the 7 AT — > n a 7rpN reaction is given 
by 

1 M N 1 f d 3 p 2 d 3 q a d 3 q p M N 1 1 



J_ f d 3 p 2 d A q a d A q 
, rel ] (2tt) 3 (2tt) 3 (2tt 



2|k| E P1 v rel J (2tt) 3 (2tt) 3 (2tt) 3 E P2 2u a 2uj p 



x (2rr)W\ Pl + k-p2-q a -qp)J2 \\ (1/2, A T |l/2, v) | 2 , (13) 



2' 

vv' 
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where k— (|k|, k), pi = (E Pl ,p 1 ), p 2 = (E P2 , p 2 ), and g 7 = (u; 7 , q 7 ) (7 = a, (3) are the four- 
momenta of the initial photon, initial nucleon, final nucleon, and outgoing pion, respectively, 
v re i is the relative velocity between the initial nucleon and photon, and u; 7 = yj m\ + q 7 and 
E p = a/M^ + p 2 . The cross section is evaluated in the 7iV center-of-mass system. The T 
matrix is in general expressed as T = A + icr ■ B, which is summed over the final nucleon 
spin states (z/) and averaged over the initial nucleon spin states (u). 

We show how to evaluate the T matrix by taking one of the diagrams as an example. Let 
us consider the process corresponding to (a) in Fig.l computed with the PS coupling. The 
T matrix is divided into two parts of the hadronic process and electromagnetic process and 
then into the particle and anti-particle intermediate states for the convenience. Thus, the 
T matrix is expressed as 



y^Ty-H {1) 1 ^#(1) 

A 1 A — n-yNN p _ U _ p p -"JVi^ 



+ n-yNN p _ 7 , p p n NNTTi^ 

where the transition matrices H are given by 

H jnn = eM(p 2 )(i 71 i7 • e - iF 2 a tv e i k l/ )u(p 2 - k) 

*X*e(S$ + itr-VV>)x, (15) 



with 



0(1) _ Je^+m£ Je^+m^ ( F 1 + F 2 k F\ - F 2 k \ 

S e -\/^M^\/^M^[ Ep2 _ k + MN + E P2 + M N ) P2 - e ' (16) 

v (i) _ / e P2 +m n / E P2 - k +M N r F\ + {Epz-k + M N + k)F 2 

E V 2M N V 2M N I ^_ fc + M N 

F X + F 2 A; F 1 -F 2 k\ 

P2 x e , (17) 

E p2 —k + iWjv E P2 + M N J y2 J ' 1 ; 



for the 7./V — > iV process, 



rr(l) -/ n 5 i>l-ia + M N 5 
^JVATtttt = M (P2 - k)^Aw7 ^ TTTQirNNl T a U{ Pl ) 

yPi la) J-yJ-N 



^X^NNisP+Ur-VP)* (18) 
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with 



_ . /E P2 _ k +M N /e p1 +m N[/ _ / q« • pi q«-(pi- k) 



1 V 2M^ V 1M N L"« ^ +Mn + ^ + Mn 



1 1 

£ P2 - fc + M ~ + M 



y(l) _ . E P 2 -k+M N E P1 +M N r I , , x 



£ p2 _ fc + AV ( Pl -g Q )2_ M 2 > 

for the iV — > 7r7riV process, 



with 



for the 7 — > NN process and 



with 



,.(2) E p k +M N I E Pl +Mjy r . / 1 1 \ , „ 

+ Qj0)i T (i T a 



E p2 _ k + M N l ( Pl -q a y-M 2 / 



(20) 



#$jv = «(P2)e(F l7 • e - iF 2 <r i, 'e i A; I/ )u(-p2 + k) 

X^-O^+^'Vg^x, (21) 



*P = (^f^ + j (k x .) • p 2 , (22) 



+ F2 '^T^-^ 2 T^» (kX£)xp2] ' (23) 



o-(2) -/ , I \ 5 • A - rfa + M 5 / \ 

^AWtttt = V{-P2 + k)c/ 7rA rjv7 T /3 ^ pi _ " )2 _ M 2 &riW7 T-aU(Pl) 

^xVrNNi-WsP + ia-V^x, (24) 



5{ 2) = 0, (25) 



(26) 



J P2-k t 1V *N KFL Ha) 

for the NN — > 7r7r process, respectively. Here u and t> are the Dirac spinors for nucleon and 
anti-nucleon, respectively and x is the two-component spinor. In this expression, we neglect 
the 0((E P + M N )~ 2 ) contributions. However, we found this approximate T matrix gives 
nearly the same result as the exact one within the present energy region. The calculations 
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for the other diagrams are performed in a similar fashion. In this sense, the T matrix is 
evaluated in a relativistic way. 

In order to estimate the relativistic effect, we calculated the cross section with the PV 
coupling for the diagrams (a)-(g) in Fig.l in both relativistic and non-relativistic ways. Here 
non-relativistic approximations mean that both the anti-particle contributions and terms of 
order (p/M^) 2 or higher in vertex operators are neglected but the denominators of the 
propagator are treated in a relativistic way. We found that the relativistic calculations are 
about 30-60% larger than the non-relativistic calculations around the iV*(i520) resonance 
energy. This indicates that the relativistic effect, which have been so far neglected in previous 
studies, should not be discarded for the non-resonant process in the 7iV — > nirN reaction. 

Now we show the calculations of total cross sections for three isospin channels: 7£> — > 
7r + 7r~p, 7p — > 7r + 7r°n and ^p — > n n°p in Fig.3. Here solid and dashed lines correspond 
to the full non-resonant calculations for the PS and PV models including the vector meson 
contributions (Fig.2), respectively. We observe that the PV results of the ir + ir~ channel are 
close to the PS results, while in the ir + ir° and 7r°7r° channels there are significant differences 
and the PS calculations are larger than the PV calculations. The similar feature can be seen 
in the one-pion photoproduction where the difference between two couplings is prominent 
in the neutral pion photoproduction but is very small in the charged pion photoproduction. 
This arises certainly from the different off-shell behavior between the PV and PS irNN 
couplings. For the one-pion photoproduction, such difference appears in the Born terms 
proportional to the anomalous magnetic moment obtained using the second term of Eq.(4). 
In fact, if F 2 is set to zero, the difference between the PS and PV calculations disappears. The 
same things happen for the two-pion photoproduction, which has been examined numerically. 

It is interesting to observe that the PS calculation in the ^p — > 7i + 7i°n channel is re- 
markably larger than the PV calculation and is roughly consistent with the size of the p 
Kroll-Ruderman term introduced in the model of Ochi et al.[12, 13], by which the large 
measured cross sections have been successfully explained. In order to reproduce the data, 
however, the small range-paramater of the pim form factor must have been used and it gave 
rise to the larger tttt scattering p-wave phase shift at low energies compared with the exper- 
imental value. We think that the p Kroll-Ruderman term with this form factor simulates a 
non-negligible background process originating from the strong coupling between the nucleon 
and I = J = 1 nn system as is inferred from strong enhancement at low energies in the 
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FIG. 3: Non-resonant contributions to total cross sections for 7p — > 7r + 7r~p, — ► 7r + 7r°n and 
■7P — > 7r°7r°p reactions. The solid and dashed lines correspond to the PS and PV calculations with 
the vector meson contributions, respectively. The meaning of remaining lines is given in the text. 

isovector rnr spectral function derived from the irir — > NN helicity amplitude [15]. Therefore, 
the p Kroll-Ruderman term is considered to represent the non-resonant process arising from 
the PS coupling effectively. 

In Fig. 3, we also show the results for the Born terms coming from the PS coupling (dotted 
lines) corresponding to the diagrams (a)-(e) in Fig.l. We find that the contribution to the 
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FIG. 4: Non-resonant contributions to the 7ttt invariant mass spectra for yp — > tt + tt p, jp — > ir + 7r n 
and 7p — ► 7T tt £> reactions at 750 MeV. The meaning of lines is the same as that of Fig.3. 



yp — > 7T + ?T~p channel is extremely large and is mostly attributed to the pion pole terms 
(diagrams (d) and (e) in Fig.l). The similar situation occurs when one calculates the nN 
isoscalar s-wave scattering length with the PS coupling. These unfavorable results can be 
improved by introducing the contact interaction of Eq.(3) as mentioned before. This effect 
is seen in the calculations with the contact terms (dash-two-dotted lines in Fig.3). Here 
the contact terms correspond to the diagrams (h)-(k) in Fig.l and consist of the isoscalar 
and isovector parts. The isoscalar term contributes to the •yp — > n + 7i^p and yp — > 7v 7i°n 
channels and the isovector term contributes to the •yp — > n + n^p and yp — > n + 7T°n channels, 
and the isoscalar term has larger coupling constant than the isovector term as understood 
from Eq.(3). The size of these terms can be seen in the calculations in Fig.3. 

For comparison, the result for the Born terms coming from the PV coupling (diagrams 
(a)-(g) in Fig.l) is also shown by dash-dotted line in Fig.3. The dash-dotted line is drawn 
only in the yp — > tt + tt~p channel but omitted in other channels since the full non-resonant 
calculations (dashed line) in other channels are almost overlapped with the dash-dotted 
line. The difference between the solid and dash-two-dotted lines or between the dashed and 
dash-dotted lines arises from the vector meson contributions of the diagrams in Fig.2. The 
contributions are relatively small. This smallness is mainly attributed to the effect of form 
factors used. In this paper, we do not pursue this effect furthermore. 
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In order to see the difference between the PS and PV models in further detail, the tttt 
invariant-mass spectra are calculated and the results at 750 MeV are shown in Fig. 4. The 
meaning of lines is the same as that of Fig. 3. Apart from the size of the distributions, the 
difference of the shape can be seen in these distributions, especially in the 'jp — > 7r°7r°p 
channel. Even in the PS calculations for the 7p — > 7T + 7T~p and 7£> — > 7r + 7r°n channels, there 
are some shifts of the peak to the higher invariant mass compared with the PV calculations. 
Here it is interesting to note that the I = J = 1 nn system relevant to the p meson can 
contribute to the jp — > n+n^p and 7p — > 7r + 7r°n reactions and on the other hand the 
/ = J = 7T7T system relevant to the a meson can contribute to the jp — > 7i + n^p and 
7P — > 7r°7r°p reactions. The correlations for the final system might influence both the 
shape and size of the distributions but in our calculations they are not taken into account. 
In this work, as a first step, we would like to demonstrate how different the PS model and 
PV model are within our present framework. In order to compare the calculations with the 
experimental data, the resonant contributions must be included. We will employ a simple 
model for the resonant processes, which will be discussed in next section. 

III. RESONANT PROCESSES 
A. isobar model 

The two-pion photoproduction in the second resonance energy region involves the res- 
onances such as p meson, A(1232) and iV*(1520) as important intermediate states. We 
treat these resonances with the isobar model where the scattering of nn or nN in the rele- 
vant channel is assumed to be described solely by the resonant state. We employ the same 
formalism used in Ref. [12, 13], which is briefly reviewed in this section. 

In this model, the 7T7T p-wave scattering t matrix in the energies from threshold to the p 
meson resonance is assumed to be written as 



t. 



(27) 



TTTT 
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Fpn-K — 2hp(n)(£ p ■ k), 

where f pn7T and q p7T7T are the pun coupling constant and range parameter, respectively and 
e p is the polarization vector and k is the relative momentum between two pions. The p self- 
energy S p7rvr is evaluated with the same way used in Ref.[12]. The parameters m° p , / p7r7r and 
q p7T7T are determined to fit the nn p-wave phase shifts as well as the mass and width of p 
meson. We take m° p = 910 MeV, / p7r7r = 7.8 and q p7T7T = 800 MeV/c, which are used in 
the calculation of the N* self-energy. In the previous paper [12, 13], q p7T7T was adjusted to fit 
the 7p — > 7T + 7r°n data instead of the nil p-wave phase shifts. In the case of Ref.[13], q p7rn 
was taken to be 200 MeV/c, which made the size of the p Kroll-Ruderman term quite large. 
This parameterization was only a way to phenomenologically reproduce the jp — > n + 7i n 
reaction cross section within the previous model. 

The n N scattering t matrix in the P 33 channel is written as 

, F-kN A-^ttJVA /OCA 

*™ - VS-A/»-EV (29) 
where \/s and M A denote the ttN center-of-mass energy and bare mass of A (1232), respec- 
tively. The vertex function for the nN — > A transition is expressed as 



Fl NA = -<v^y^5^A(p)(St • P ), (30) 

where p is the three-momentum in the nN center of mass system and p is its unit vector. 
St is the spin transition operator from 1/2 to 3/2 and g^NA is given by[18] 

9nNA{P) = == === o ' ( 31 ) 

y/2(m n + M N ) m n \Q 2 A +p 2 J 
where is the coupling constant and Qa is the range parameter. The A self-energy T,^ N is 
evaluated using the vertex function F^ NA . The parameters M A , and Qa were adjusted 
to fit the experimental P33 scattering amplitude [18]. 

The nN scattering t matrix in the D 13 channel has the same form with the the above P 33 
amplitude tp 33 as follows: 



4 

(32) 



F nNN *F nNN , 
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where M%« denotes the bare mass of N*. The vertex function for the irN — > N* transition 
is written as[19] 

Fl NN , = -.W/'j y (^)V<p/.w>' (SOU . y 2 (p)) , (33) 

V ^Wjv y/2(m n + Mat) \PttNN* J 

where f^NN* is the tcNN* coupling constant and p w nn* is the irNN* range parameter. 
is the second-rank spin transition operator from 1/2 to 3/2, which is defined by 

= y| [ S t x CT ] ( 2 ) . (34) 

The iV*(1520) resonance can decay into both the 7riV and irirN channels. The tcttN decay 
occurs through three dominant modes, i.e., (nA) s _ wave , (irA) d _ wave and pN . These branch- 
ing fractions are known to be comparable. Thus the total N* self-energy (E to tai) is expressed 
as 

^total — ^VTV + + ^VA + ^pN, (35) 

where S^, E* A , E^ A and S pA r are due to the coupling to the irN, s-wave 7rA, d-wave it A 
and pN channels, respectively. is evaluated from the vertex function Fl NN , and the 

other components of the self-energy are obtained from the following vertex functions: 



KIn*(p) = ^(2^) 3/ %/ 2(m ^ M ^ /^e-(^0V 00 (p), (36) 



for the N* -> (ttA) 



s—wavei 



F^lp) = -,2^/y_^ A „. (-/-) .-<^ ( S - ■ V' 2 (p)) , ,37) 
for the N* -»■ (7rA) d _ wa „ e and 



F p U, = (^) 3/2 /Swe-W^*) 2 (St • e p ) F 00 (p), (38) 

for the iV* — > pA 7 " . Here f^AN* , f^AN* and fpNN* are the s-wave, d-wave ttAN* and 
pNN* coupling constants, respectively. p* AA r* and p p nn* are the 7rAA^* and pNN* range 
parameters. is the second-rank spin transition operator from 3/2 to 3/2 defined in 

Ref.[12]. E^ A and S pA r contain the effect of the decay process A — > ttN or p — > tttt and 
their explicit forms are given in Ref.[12]. 
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In this paper, we simply choose 400 MeV/c for the range parameters (p n NN*, P s w an*i 
P^atv*; PpNN*) which reproduces the nucleon size in quark models[19]. The coupling con- 
stants (fnNN*, I-kAn* ' fnAN* i fpNN* ) an d the bare mass (M^r*) are adjusted to fit the N* 
resonance energy, its width and the branching ratios at the resonance energy We use 1520 
MeV as the resonance energy and 120 MeV as the width, respectively. We take a fraction 
of 60% for the decay into irN, 8% into s-wave ttA, 12% into d-wave irA and 20% decay 
into the pN channel, respectively [20], which are slightly different from the values used in 
the previous model[12, 13]. We note that the parameters cannot be uniquely fixed due to 
the limited experimental information and their uncertainties. The parameter set used in 
this paper are given in Table 1. The signs of the coupling constants are the same as those 
in Ref. [13]. With this parameterization, the nN D 13 scattering amplitudes calculated agree 
with the data around the resonance energy but are deviated as the pion energy is away from 
it. The range parameters are necessary to be varied in order to get good agreement with 
the data over a wide range of the pion energies, but they are not uniquely determined in 
the present case and the fitted parameter sets include a very small range parameter which 
is hardly acceptable from a physical point of view. Furthermore, the background nN inter- 
action may be needed and thus one must go beyond the present framework of the isobar 
model. We think that the above description for N* is sufficient for the present purpose, 
namely, to examine the difference between the PS and PV couplings. 

TABLE I: Parameter set for N* used in this paper 



fnNN* 


fwAN* 


JwAN* 


fpNN* 


M^(MeV) 


1.147 


0.398 


1.435 


0.942 


1709. 



Now we consider the resonant couplings by photon. The 'jNA and '-/NN* coupling 
constants can be determined by using the multipole amplitudes in the relevant channel 
for the 7iV — > nN reaction. The multipole amplitudes for the 7iV — > nN reaction have 
non-negligible background contributions and thus are generally expressed as the sum of the 
background and resonant terms, i.e., 

t1 n = T f + T f. (39) 
Like the nN elastic scattering amplitudes, the resonant term T^ N is given by the isobar 
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model as 

Tl N _ F 7tNA F} /NA 

for the A resonance and 
for the N* resonance. 

The A resonance can contribute to both M 1+ (3/2) and Ei + (3/2) multipole amplitudes. 
Since the magnitude of the E 1+ (3/2) multipole is small compared with M 1+ (3/2), the E2 
>yNA coupling is neglected. The 7^ — > A vertex function for the M 1+ (3/2) channel is 
written as 



Fj fNA = -ig Ml+ &-kxe, (42) 

where Qm 1+ and e are the Ml 7iVA coupling constant and photon polarization vector, 
respectively, and k denotes the initial photon momentum. We use <?Mi+ = 0.1991 (in natural 
unit), which is obtained from the resonant coupling given by the Particle Data Group [20]. 

For the N* resonance, we use the helicity amplitudes instead of the electric and magnetic 
multipoles. The •yNN* vertex has two independent helicity couplings: the helicity 1/2 and 
3/2 couplings. For the helicity 1/2 transition, F^ NN * is written as 

F' ! nn* = -m/2 (S f • k) (<r-k x e), (43) 

where gi/2 is the helicity 1/2 coupling constant and k denotes the unit vector of the initial 
photon momentum. For the helicity 3/2 transition, F^ NN „ is written as 

F%%* = 9 V 2 { ■ e) + \ (St • k) (<r k x e) } , (44) 

where is the helicity 3/2 coupling constant. For the proton target, the helicity 1/2 
amplitude is small compared with the helicity 3/2 amplitude and so the helicity 1/2 coupling 
is neglected. For the neutron target, both couplings are taken into account. We use gi/ 2 = 
and g 3 / 2 = 0.1612 for the proton target and gi/ 2 = —0.0496 and g 3 / 2 = —0.135 for the 
neutron target, respectively, which are obtained from the resonant couplings given by the 
Particle Data Group [20]. 
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Generally, the 'jNA and 'jNN* coupling constants ( gn 1+ , 9i/2 , <?3/2 ) include both 
bare couplings to the resonances and vertex corrections due to the interference with the 
background processes and thus are complex and energy dependent in nature. One way 
to determine the coupling constants is to extract them from the experimental multipole 
amplitudes by assuming an appropriate non- resonant background term T^ N . The other way 
is to use the resonance couplings given by the Particle Data Group [20] which correspond 
to the bare couplings. In most of previous models[10, 11, 14], real coupling constants 
obtained by the latter method have been used to calculate the cross sections for the two- 
pion photoproduction. Since the imaginary values are small for these resonances as predicted 
in the phenomenological calculations [16], we also employed the latter method for simplicity. 



B. resonant amplitude of 7^ — > ttttN 

We use the same approach with the model of Ref. [12, 13] for the resonance production 
processes and make a brief review about it in this subsection. The resonant Tr matrix for 
the two-pion photoproduction is expressed as 

Tr = T/S.KR + TapP + 7yv*7rA + ^N"*7rA + Tn* p N + T^nex- (45) 

The resonant Tr matrix consists of six amplitudes: the A Kroll-Ruderman term (T^kr), 
A pion-pole term (T APP ), N* excitation terms (T^?i A and T N * pN ) and the 7rA production 
term accompanied by nucleon exchange (T& nex ). These diagrams are shown in Fig. 5 (a)- 
(e), respectively. The strength of these resonant processes depends strongly on the isospin 
channel[10, 12]. The terms T^kr and T^pp dominate for the -k + tt~ production and becomes 
small due to the isospin factor for the 7r + 7r (7r _ 7r ) production and then are prohibited for 
the 7T 7r° production. Although the N* excitation terms have only weak strength, they 
contribute to all isospin channel and their interference with the A Kroll-Ruderman term has 
significant effects to the cross section. The other processes arising from the requirement of 
the gauge invariance are neglected because they are known to be small. 
The A Kroll-Ruderman term is written as 



FnNAF ai 

v/s - uj n - E A - J: { l N) (q, y/s) 



Takk = — ™ "Sim. ^ ( 46 ) 
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(a) (b) 




(d) (e) 

FIG. 5: Diagrams of resonant processes. (a) the A Kroll-Ruderman term (Takr), (b) A pion-pole 
term (Tapp), (c) N* excitation terms (?^f| A ), (d) N* excitation terms (Tn* p n) and (e) the ttA 
production term accompanied by nucleon exchange (T& nex ). 

where uj w = y/rn^ + q 2 and = a/ (M^) 2 + <? 2 - ^a^C?) \A) ^ s the self-energy of A moving 
with the momentum q that arises from the coupling to the ttN channel. Its expression is 
given in Ref. [19]. The 'jNnA contact term operator F AKR is obtained from the strong 
7riVA vertex function so as to satisfy the gauge invariance, the detail of which is given in 
Ref. [12]. Instead of using the effective Lagrangian[10], we employ the vertex function with 
a form factor. The N — > 7rA transition operator is assumed to have the same form as the 
A — > tvN vertex function. Since the range parameter Q&(N — > 7rA) is not necessarily the 
same as the parameter Qa in Eq.(31) and thus is unknown, we treat it as a free parameter 
and vary it to fit the — > n + n^p cross section. The pion-pole term T^pp is obtained by 
replacing -F AKR in Eq.(46) with the 'jNirA pion-pole vertex function F APP whose detailed 
expression is given in Ref [12]. 
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The 7iV — >■ N* — > nnN transition takes place the following two possible processes: 
N* — > 7i A (s-wave or d-wave 7rA state) and N* — > piV. They are described by and 
Tn* p n, respectively. Using the isobar model mentioned in the previous subsection, the T 
matrix elements of T^fl A and T N * pN are written as 

TP p s Wt pt 

j,s(d) _ r ^NA^ n AN* r ^yNN* ,^„x 

rp _ FpnnF p NN"F^ NNt 

N * pN ~ 2u p (^s-u p -E qp -i: p Uq P ^)){^-M% t -^y 1 ] 

respectively. Here, q p is the momentum of p meson and S p7r7r (g p , y/s) is the self-energy of p 
meson moving with the momentum q p . The p meson is described by the isobar model (see 
Eq.(27)). 

The 71 A production term accompanied by nucleon exchange is written as 

TP pt pt 

rp^ _ r ^NA r -yNA r nNN 

( y/s - ui n - E A - S^ rAr) (q, y/s) ) (E k -u n - E k+q ) ' 
where the ttNN vertex function F^ NN is assumed to have a usual non-relativistic form. 
Unlike the other resonant processes, the intermediate particle, i.e., nucleon, is far off-shell. 
Since only the on-shell resonant Mi + (3/2) multipole amplitude is known, some prescription 
is needed to include the off-shell effect in the above T matrix (Eq.(49)). We take the modified 
pole approximation [21]where the angular part of i ? j iVA is evaluated at the center of mass 
system of 7iV and its magnitude is evaluated at the total energy of the final ixN state. The 
contribution of T& nex to the cross sections is expected to be small due to the off-shell effect. 
However we include this amplitude in our calculations, since it becomes non-negligible for 
the 7p — > 7r n°p reaction due to large coupling constants in Mi + (3/2) channel. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we present our numerical results of total cross sections, invariant mass 
spectra and helicity-dependent cross sections obtained by using the model introduced in 
the previous sections and compare them with the experimental data. In our model, the T 
matrix for the two-pion photoproduction is written as 

T = T NR + T R , (50) 
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where T^r is the non-resonant T matrix given in section II and Tr is the resonant T matrix 
(Eq.(45)) in section III. For T NR) we consider two kinds of models: PS and PV models 
with vector meson contributions. Correspondingly there are two kinds of full T matrices for 
Eq.(50) which are expressed as Tp$ and Tpy, respectively. To demonstrate the difference of 
nNN couplings, we always compare the calculations by T PS with those by T PV which are 
shown in figures below. 

In our model, there is a free parameter: the range parameter Qa(N — > nA) of the 
N — > nA vertex function appearing in the A Kroll-Ruderman term. At present we do 
not know how to determine it by using other reactions than the two-pion photoproduction. 
We use Qa(N — > 7rA) =430 MeV/c determined so as to reproduce the 7£> — > n + n~p cross 
sections. Here the N — > 7rA vertex function is assumed to have the same form with the 
A — > nN vertex function used in this paper, since the A propagator used in Eq.(46) is 
calculated with the same A — > nN vertex function. In this sense, we think, it is not 
appropriate to use the monopole type of form factor employed by the authors of Refs.[10, 14] 
in our approach. Although the value of Qa(A — > nN) is taken to be 358 MeV/c given by 
Betz and Lee[18], the value of Qa(N — > nA) does not need to be the same, since the form 
factors for A — > nN and iV — > nA are functions depending on the relative momentum of 
7i N and ttA systems, respectively. At a given relative momentum, in fact, the square of the 
pion four-monentum evaluated in the nN center of mass system is larger than that in the 
7rA center of mass system due to the mass difference between N and A. 

First, we show the results of the total cross sections of five isospin channels, i.e., ^p — > 
n + n~p, 7p — > n + n°n, ^p — > 7r°7r°p, — > n + n~n and — > n~n°p, in Fig. 6. Here the 
calculations with Tp$ and Tpy are presented with the solid and dashed lines, respectively. 
We find that the PS calculations are in good agreement with the data in all channels and on 
the other hand the PV calculations underestimate the data except for n + n~ channels. The 
discrepancy between two calculations is clearly due to the difference of the non-resonant 
processes. Our results of the 7p — > n + n°n and 7^ — > n°n°p reactions for the PV model 
are not consistent with those by the model of Nacher et al.[14] where the non-resonant 
amplitudes are constructed using the non-relativistic PV nNN coupling. Their results are 
more close to the experimental data compared with our result for the PV model. We think 
that this inconsistency comes mainly from the way of the calculation of the diagram (e) in 
Fig. 5, where the nucleon in the intermediate state is far off-shell. We have evaluated this 
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FIG. 6: Total cross sections of two-pion photoproduction on proton and neutron in various isospin 
channels. The solid and dashed lines correspond to the PS and PV calculations. Data are from 
Refs.[l] (triangle down), [2] (open square) , [4] (black diamond), [5] (star), [6] (open circle), [7] (triangle 
up and black square), [8] (open diamond) and [9] (black circle) . The data (black circle) of 771 — > 
Tr~ir p correspond to the cross section over the DAPHNE acceptance. 

diagram with the modified pole approximation and on the other hand they adopted simply 
the 7 TV A vertex function without cut-off factor which has a linear dependence of the photon 
momentum. Consequently, the 7./VA coupling becomes very large in the energy region of 
the iV*(1520) resonance and leads to the large cross section. 

Secondly, we calculated the invariant mass spectra of three isospin channels, i.e., 7^ — > 
7i + 7i^p, 7£> — > 7i + 7i n and 7p — > n°7i°p. The invariant mass spectra for the tc + tt~ system 
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FIG. 7: Invariant mass spectra of irn for — > 7r + 7r~p at 650, 750 and 850 MeV. The solid 
and dashed lines correspond to the PS and PV calculations, respectively. Data are normalized 
appropriately and from Ref.[l]. 



in the ^yp — > 7i + 7r~p channel are shown in Fig. 7. The calculations are performed with 
three photon energies, i.e., 650, 750 and 850 MeV. In this case, the values of the data are 
plotted arbitrarily and the calculations are normalized so as to fit the peak value of the 
experimental distributions. For the shape, two calculations (PS and PV calculations) are 
almost equivalent each other and are in good agreement with the data. 

For the •yp — > 7r + 7r°n channel, the invariant mass spectra for the 7r + 7r°, n°n and n + n 
systems are calculated at four bins of photon energies, i.e., 650-700 MeV, 700-740 MeV, 
740-780MeV and 780-820 MeV and are shown with the data in Fig.8. We find that the 
PS calculations agree well with all data but the PV calculations show some discrepancy 
about the shape of the 7r + 7r° invariant mass spectra in addition to the magnitude. The 7r + 7r° 
invariant mass spectra have a peak shifted to the higher nir invariant mass and the peak 
position of the n°n and n + n invariant mass spectra corresponds to the mass of A(1232). In 
order to see which process makes such behavior, the contributions of some components of the 
resonant processes are shown in Fig. 9. Here the dash-dotted and dashed lines correspond to 
the calculations with T^* p n and Takr + Tapp + T& nex + T^* nA + T^* nA , respectively. For 
reference, the non-resonant PS calculations (solid lines) are also plotted in figures. Clearly, 
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FIG. 8: Invariant mass spectra of tt + tt°, ir n and 7r + n for 7^ — > 7r + 7r°n at four bins of incident 
photon energy. The solid and dashed lines correspond to the PS and PV calculations, respectively. 
Data are from Ref.[6] 
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FIG. 9: Contributions of resonant processes in invariant mass spectra of 7r + 7r°, ir°n and 7r + n 
systems for — > 7r + 7r°n at 750 MeV. The dash-dotted and dashed lines are the calculations with 
T N * pN and Takr + Tapp + T& nex + T^» 7rA + 7$,^ (see text), respectively. The non-resonant PS 
calculation (solid line) is also plotted. 

the peak shift for the ir + ir° invariant mass distribution is due to the contribution of N* — > pN 
process and the peak position of the n n and n + n invariant mass distribution is directly 
related to the A(1232) production in the intermediate state. 

For the 'jp — > TT ir°p channel, the invariant mass spectra for the 7r°7r° and ir°p systems are 
calculated at the same bins of photon energies with the 7p — > n + 7i n channel. The results 
are shown with the data in Fig. 10. In this channel, the A Kroll-Ruderman term (Fig. 5 (a)) 
and pion pole terms (Fig.5(b)) do not contribute to the cross section and so the magnitude of 
the cross section is rather small compared with other channels. Furthermore the I = J = 1 
7T7T system such as the p meson is not produced because of isospin conservation. Therefore 
only the processes of T / \ nex ,T^* nA and T^-* 7rA among the resonant processes contribute to the 
cross section in our model. We note that the production of the I = J = rnr system such 
as the a meson could take place in this channel. From the comparison with the data, one 
finds that both PS and PV calculations can almost reproduce the data of the ir°p invariant 
mass spectra, which has a peak at the A(1232) mass. For the 7r°7r° invariant mass spectra, 



27 




FIG. 10: Invariant mass spectra of tt it and ir p for 7p — > 7r 7r p at four bins of incident photon 
energy. The solid and dashed lines correspond to the PS and PV calculations, respectively. Data 
are from Ref.[8] 
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there are some discrepancies, especially, between PS calculations and the data. The PV 
calculations are almost consistent with the data except for the bin of 780-820 MeV, where 
there are two bumps in the distribution. On the other hand, the PS calculations show 
broader distributions than the data, arising from the non-resonant processes as shown in 
Fig. 4. In our calculations, the final state interactions for 7r°7r° and 7i°p are not taken into 
account. As pointed out in Ref. [22], the final state interaction for 7r°7r° is expected to be 
important because of strong correlation between two pions in / = J = channel and the 
influence of the a meson. Such effect might be one of the possibilities to improve the PS 
model and so it will be worthwhile to further study this channel by taking into account such 
effect. 

Finally the calculations of the helicity-dependent cross sections 03/2 and 0i/ 2 are shown 
with the data in Fig. 11. Here we display the cross sections in three isospin channels: the 
7P — > it + it~p cross section with the three charged particles in the DAPHNE acceptances [23, 
24], the 7p — > TT + -K°n cross section with the ir + in the DAPHNE acceptance [23, 24] and the 
7P — > 7i 7i°p cross section without kinematical limits for outgoing particles. The quantity 
°3/2 (ci/2) is defined as the cross section for the absorption of a polarized photon by a 
polarized target proton in the helicity 3/2 (1/2) channel. The solid and dashed lines denote 
the calculations with the PS model and PV model, respectively. It should be noted that the 
electromagnetic coupling of the iV*(1520) to proton is known to be dominated by the helicity 
3/2 state[20]. This is why the peak of the resonance is seen in the data for 0-3/2 but not for 
01/2. Therefore the cross section 01/2 is sensitive to the reaction mechanism not related to 
the iV*(1520) resonance. For the ^p — > -K + n~p reaction, one finds that the PS calculation 
is almost equivalent with the PV calculation like the total cross sections. Although the 
calculations agree qualitatively with the data, the cross section 3 / 2 is underestimated and 
the cross section 0^2 is overestimated. The detailed description of the reaction mechanism 
is still unsatisfactory for this channel and is needed to pursue what is missing in our model. 
For the jp — > 7i + n°n reaction, the PS model can almost explain two helicity-dependent cross 
sections simultaneously except for the cross section ay 2 around 750 MeV. On the other hand, 
the PV model underestimates the cross section 3 / 2 and also 01/2 slightly. In this model, 
non- resonant processes in the helicity 3/2 channel are not strong enough to explain the data. 
This result is not consistent with that obtained from the model by Nacher et al.[23]. This 
is due to the same reason mentioned in the discussions of the total cross sections. For the 
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FIG. 11: Helicity-dependent cross sections a 3 / 2 and a x / 2 for 7^ — > 7r + 7r~p, 7^1 — ► 7r + 7r°n and 
7P — > 7r°7r°p- The solid and dashed lines correspond to PS and PV calculations, respectively. The 
upper two lines and black circles correspond to the cross section <t 3/ / 2 and the other two lines and 
triangles correspond to the cross section ai/ 2 in each figure. In the above two figures of — > ir + 7r~p 
and 7p — > 7r + 7r°n , the cross sections over the DAPHNE acceptance are shown. Data are from 
Ref.[23, 24] 

7P — > 7r°7r°j9 reaction, we observe that there is a large difference between the PS and PV 
models in the cross section a 1 / 2 . The experimental data, if exist, could provide an important 
information on the reaction mechanism. 
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V. CONCLUDING REMARKS 



Several theoretical studies on the two-pion photoproduction have been performed in 
the past but none of them have succeeded to reproduce the data in all isospin channels 
simultaneously. Particularly, unexpectedly large cross sections of the ^p — > n + n°n and 
7n — > 7T~7T°p reactions were found not to be explained with the usual reaction mechanism 
and thus the presence of a new reaction mechanism in these channels is suggested. Ochi 
et al.[12]introduced the p Kroll-Ruderman term, which influences only on the above isospin 
channels, as a new reaction mechanism. However, a soft pnix form factor was needed to 
reproduce the large cross sections and the isobar model with such a form factor failed to 
explain the mi scattering at low energies. Thus, the p Kroll-Ruderman term was inferred to 
represent the effect of the background process followed by the production of the I = J = 1 
nn system rather than the p meson itself. 

In this paper we have discussed the effect of the non-resonant processes arising from 
the PS and PV ttNN couplings in the two-pion photoproduction in order to pursue an 
alternative reaction mechanism. The non-resonant amplitudes can be obtained by attaching 
an external photon line to the diagram of iV — > irirN calculated with the two irNN couplings. 
In order that the model with the PS irNN coupling is equivalent to the model with the PV 
iiNN coupling in the strong interaction process N — > ttttN, the effective contact interaction 
(Eq.(3)) is added to the former model. Using these models, we examined the effect of the two 
couplings on various isospin channels numerically. For the jp — > 7i + n°n and — > 7v^iv p 
reactions, we found that the cross sections calculated with the PS model were larger than 
those with the PV model and their magnitude was almost consistent with the p Kroll- 
Ruderman used in Ref.[12]. Consequently, the non-resonant process described by the PS 
nNN coupling can be regarded as a candidate for the new reaction mechanism in place of 
the p Kroll-Ruderman term. 

The importance of the PS nNN coupling in the two-pion photoproduction has not been 
so far noticed. In fact the non-resonant process has been always described by the PV 
ttNN coupling and has been found to have only a minor contribution to the two-pion 
photoproduction. The PV ttNN coupling is more favored than the PS nNN coupling for 
the 'jN — > ttN reaction and the nN scattering at low energies. Therefore, the PV coupling 
has been also used for the ^N — > nnN reaction as a matter of course. Recently, Drechsel 
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et al. [16] pointed out in the study of the 'jN — > irN reaction that the PS coupling was 
needed to describe the data with the increase of the incident photon energy. This implies 
that the PS coupling is preferable to the PV coupling at larger off-shell nucleon momenta. 
Accordingly the PS coupling is expected to be important for the two-pion photoproduction, 
since the far off-shell nucleon is involved in the intermediate state. 

In order to demonstrate the importance of the PS coupling and compare our theory 
with the data, we constructed two types of models (PS and PV models) by adding the 
vector meson contributions for the non-resonant processes (Fig.2) as well as the A(1232) and 
iV*(1520) contributions for the resonant processes (Fig.5) to the non-resonant contributions 
coming from the 7iNN couplings (Fig.l). Using these models, we calculated the total cross 
sections, invariant mass spectra and helicity-dependent cross sections for various isospin 
channels. Generally, these observables are successfully described by the PS model compared 
with the PV model except some details. Particularly the calculations of total cross sections 
and invariant mass spectra for the jp — > 7i + n n and 7n — > 7r^n p channels by the PS model 
are in good agreement with the data and on the other hand those by the PV model are 
largely underestimated. There are some discrepancies between the calculations with the PS 
model and the data in the invariant mass spectra for the 7p — > 7r°7r°p channel. Because the 
magnitude of the cross section in this channel is very small, the higher order processes might 
emerge in the detailed structure. Therefore, we think, it may be necessary to investigate 
the mechanism originating with the a meson and the final state interaction which are not 
taken into account in our model. Although there are still unsatisfactory points, we conclude 
that the PS coupling certainly plays an important role in the two-pion photoproduction, 
especially, the 7£> — > 7r + 7r°n and — > 7T~7i°p channels. 

Now we would like to make some remarks about our model. First, the difference between 
the PS and PV couplings prominently appears in the neutral pion productions such as 
7P — > ir°p and ^p — > n + 7r n reactions. This is understood from the fact that it stems 
from the anomalous magnetic moment term in the "fNN coupling (Eq.(4)) [25]. In fact the 
difference disappears if F 2 is set to zero. In our model, F 2 is taken to be the on-shell value. 
It would be interesting to find out the influence of the off-shell effects in the 'jNN vertex 
for the two-pion photoproduction and furthermore to know how the PS coupling at high 
photon energies is connected to the PV coupling at low energies in terms of the off-shellness 
of the intermediate nucleon. We note that the off-shell structure of the >yNN vertex appears 
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as the modification of the anomalous magnetic moment term[26]. 

Secondly, the resonant processes in the present model have been treated in a naive way. 
The strong coupling constants for the A r *(1520) resonance are determined from its total 
width and the branching ratios given in Particle Data Group [20] and their form factor ranges 
are assumed to be consistent with the nucleon size in quark models. The electromagnetic 
couplings of the iV*(i520) to nucleon are taken from the PDG and their values are real 
number. In order to investigate the role of the resonance in detail, the theory must be 
refined so as to describe three reactions, i.e., irN — > irN, >yN — > ttN and 7^ — > tttcN 
, simultaneously. For instance, the effective coupling constant of 7 iViV* (1520) becomes 
complex due to the interference with the background process, although its effect is not 
expected to be large because of the small imaginary number [16]. 

Finally, we have not included the A(1700) resonant process in our model. Nacher et 
al.[14]have pointed out that the A(1700) resonance has a significant contribution on the 
7iV — > -k-kN reaction cross section due to the strong interference between the A(1700) 
resonant process and the A Kroll-Ruderman process although the contribution of only the 
A(1700) resonant process is almost negligible. They used the real 7iVA(1700) coupling 
constant given in PDG like the iV*(1520) resonance in our model. However, it was predicted 
in the phenomenological analysis of the — > nN reaction[16] that its imaginary part 
was comparable to its real part, which is quite different from the case of the ^NN* (1520) 
coupling. This complex coupling constant, if used, should largely influence the results of 
Ref.[14] about the interference effect. Further investigation is needed to draw a definite 
conclusion on the A(1700) resonant process. 
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